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Inspired by the present experimental status of charged charmonium-like states Zjf, the tetraquark 
states [cu] [cd] are systematically studied in a color flux-tube model with a multi-body confinement 
potential. The investigation indicates that charged charmonium-like states Z+(3900) or Zjf(3885), 
Z+(3930), Z+(4020) or Z+(4025), Z+ (4050), Z^~(4250), and Z+(4200) can be uniformly described 
as tetraquark states [cubed] with the quantum numbers n 2S+1 Lj and J p of l 3 Si and 1 + , 2 3 <Si 
and 1 + , 1 5 iS I 2 and 2 + , l 3 Pi and 1~, l s Di and 1 + , and l 3 Pi and 1 + , respectively. The predicted 
lowest charged tetraquark state [cu][cd] with 0 + and 1 1 .So has a energy of 3780 ± 10 MeV in the 
model. The tetraquark states are compact three-dimensional spatial configurations similar to a 
rugby ball, the higher orbital angular momentum L between the diquark [cu] and antidiquark [cd], 
the more prolate of the states. The multibody color flux-tube, a collective degree of freedom, plays 
an important role in the formation of those charge tetraquark states. However, the two heavier 
charged states Z p (4430) and Z p (4475) can not be explained as tetraquark states [cu] [cd] in this 
model approach. 

PACS numbers: 14.20.Pt, 12.40.-y 


I. INTRODUCTION 


Quantum chromodynamics (QCD) is nowadays widely 
accepted as the fundamental theory to describe hadrons 
and their interactions. Conventional hadrons are com¬ 
posed of either a valence quark q and an antiquark q 
(mesons) or three valence quarks (baryons) on top of the 
sea of qq pairs and gluons. One of the long standing chal¬ 
lenges in hadron physics is to establish and classify gen¬ 
uine multiquark states other than conventional hadrons 
because multiquark states may contain more informa¬ 
tion about the low-energy QCD than that of conven¬ 
tional hadrons. In the past several years, a charged char- 
moniumlike Z+ family, Z+(4430), Z+( 4050), Z^"(4250), 
Z+(3900), Z+(3885), Z+(3930), Z+(4020), Z+(4025), 
Z+(4475) and Z+(4200), has been successively observed 
by experimental Collaborations M- Obviously, those 
charged charmoniumlike states go beyond conventional 
cc-meson picture and prefer to tetraquark systems ccud 
due to carrying one charge, which provides a good place 
for testing various phenomenological research methods of 
hadron physics. On the theoretical side, a large amount 
of work has been devoted to describe the internal struc¬ 
ture of these charged states, which has been related 
to meson-meson molecules (TT], [l2j, diquark-antidiquark 
states [l3| . lradrocharmonium or BormOppenheimer 
tetraquarks fl3l. coupled channel cusps 0 , and kine¬ 
matic effects [lq . 

A systematic understanding of the internal structure of 
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these charged states may not only contribute to provide 
new insights to the strong interaction dynamics of mul¬ 
tiquark systems and low-energy QCD but also provide 
important information on future experimental search for 
the missing higher orbital excitations in the Z+ family. 
This is the goal of the present work. In our approach, 
a phenomenological model, color flux-tube model with 
a multi-body confinement potential instead of a two- 
body one in traditional quark model, is employed to ex¬ 
plore the properties of excited charged tetraquark states 
ccud systematically. The model has been successfully 
applied to the ground states of charged tetraquark states 
[Qq] \Q'q' 1 ( Q , Q' = c,b and q, q' = u, d, s ) in our previous 
work [171 ] ■ 

This work is organized as follows: the color flux-tube 
model and the model parameters are given in Sec. II. 
The numerical results and discussions of the charged 
tetraquark states are presented in Sec. III. A brief sum¬ 
mary is given in the last section. 


II. COLOR FLUX-TUBE MODEL AND 
PARAMETERS 


The details of the color flux-tube model basing on tra¬ 
ditional quark models and lattice QCD picture can be 
found in our previous work [l8j , the prominent charac¬ 
teristics of the model are just presented here. The model 
Hamiltonian for the state [cu][ai] is given as follows, 
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T c is the center-of-mass kinetic energy of the state, p, 
and m,; are the momentum and mass of the *-th quark 
(antiquark), respectively. The codes of the quarks (an¬ 
tiquarks) c and u (c and d) are assumed to be 1 and 2 
(3 and 4), respectively, their positions are denoted as iq 
and F 2 (r 3 and r 4 ). 

The quadratic confinement potential, which is believed 
to be flavor independent, of the tetraquark state with a 
diquark-antidiquark structure has the following form, 

V G = K [(ri - y i2 ) 2 + (r 2 - y i2 ) 2 + (r 3 - y 34 ) 2 
+ (r 4 - y 3 4 ) 2 + K d ( yi 2 - y 3 4 ) 2 ] , ( 2 ) 

The positions y 42 and y 34 are the junctions of two Y- 
shaped color flux-tube structures. The parameter K is 
the stiffness of a three-dimension flux-tube, KdK is other 
compound color flux-tube stiffness. The relative stiffness 
parameter Kd of the compound flux-tube is (l9i | 
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where Cd is the eigenvalue of the Casimir operator as¬ 
sociated with the SU{ 3) color representation d at either 
end of the color flux-tube, such as C 3 = |, Cq = and 
Cg = 3. 

The minimum of the confinement potential V G in can 
be obtained by taking the variation of V c with respect 
to yi 2 and y 34 , and it can be expressed as 


V G in = /v (R 2 + R 2 + T ^-R 2 ; , 


( 4 ) 


The canonical coordinates R, have the following forms, 


R 1 = ^=( r l- r 2 ), R 2 = ^( r 3 - r 4), 

R 3 = 4=(ri +r 2 - r 3 - r 4 ), 
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The use of V G in can be understood here that the gluon 
field readjusts immediately to its minimal configuration. 
It is worth emphasizing that the confinement V G in is a 
multi-body interaction in a multiquark state rather than 
the sum of many pairwise confinement interactions, 
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in Isgur-Karl quark model and chiral quark model with 
n = 1 or 2 . 

The central parts of one-boson-exchange V B and a- 
meson exchange V? only occur between u and d, and that 
of one-gluon-excliange V G is universal. V, B , V? and V G 


take their standard forms and are listed in the following, 
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Where x stands for n, K and 77 , Y{x) = e~ x /x. The 
symbols F, A and u are the flavor 517(3), color SU(3) 
Gell-Mann and spin 517(2) Pauli matrices, respectively. 
Op is the mixing angle between 771 and ijs to give the 
physical r\ meson. g^/Yn is the chiral coupling constant. 
a s is the running strong coupling constant and takes the 
following form [ 201 ] . 
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where jj,ij is the reduced mass of two interacting parti¬ 
cles qi (or qi) and qj (or qj). Aq, ao and /io are model 
parameters. The function d(r.y) in V G should be regu¬ 
larized I 22 L 
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where ro(^ij) = fo/iiij, fg is a model parameter. 

The diquark [cu] and antidiquark [cd] can be consid¬ 
ered as compound objects Q and Q with no internal 
orbital excitation, and the angular excitation L are as¬ 
sumed to occur only between Q and Q in the present work 
and the parity of the state [ cu] [cd] is therefore simply re¬ 
lated to L as P = (—1) L . In this way, the state [cu][cd] 
has lower energy than that of the states with additional 
internal orbital excitation in Q and Q. In order to facil¬ 
itate numerical calculations, the spin-orbit interactions 
are approximately assumed to just take place between 
compound objects Q and Q , which is consistence with 
the work a The related interactions can be presented 
as follows 
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where M 12 = M 34 = m c + m u , d , G(x ) = Y{x){± + ^), 
and S stands for the total spin angular momentum of the 
tetraquark state [ cu] [cd]. 

The model parameters are determined as follows. The 
mass parameters m n , nriK and m v in the interaction V p 
are taken their experimental values, namely, m v = 0.7 
fm -1 , rriK = 2.51 fm -1 and m v = 2.77 fm” 1 . The cutoff 
parameters take the values, A^ = A a = 4.20 fm -1 and 
A v = A k = 5.20 fm -1 , the mixing angle dp = —15° [20| . 
The mass parameter m a in the interaction V? is deter¬ 
mined through the PC AC relation m 2 ~ m 2 +4m 2 d [23| . 
m u ^ d = 280 MeV and m a = 2.92 fm -1 . The chiral cou¬ 
pling constant g c h is determined from the 7 tNN coupling 
constant through 
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The other adjustable parameters and their errors are de¬ 
termined by fitting the masses of the ground states of 
mesons using Minuit program, which are shown in Table 
I. The mass spectrum of the ground states of mesons, 
which is listed in Tale II, can be obtained by solving the 
two-body Schrodinger equation 

{H 2 ~ E 2 )$ij eson = 0. (17) 


The mass error of mesons AE 2 introduced by the pa¬ 
rameter uncertainty A Xi can be calculated by the formula 
of error propagation, 
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where Xi and Axt represent the *-th adjustable parameter 
and it’s error, respectively, which are listed in Table I. 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 


Within the framework of the diquark-antidiquark con¬ 
figuration, the wave function of the state [cu] [cd] can be 
written as a sum of the following direct products of color 
Xci isospin r/j, spin rj s and spatial <f> terms, 
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In which r, R and X are relative spatial coordinates, 
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TABLE I: Adjustable model parameters, (units: m s , m c , mi, 
po, Aq, MeV; K , MeV-fm -2 ; ro, MeV-fm; qq, dimensionless) 


Parameters 

Xi 

A xi 

Parameters 

Xi 

A xi 

m s 


511.78 

0.228 

ao 


4.554 

0.018 

m c 


1601.7 

0.441 

Ao 


9.173 

0.175 

mb 


4936.2 

0.451 

go 


0.0004 

0.540 

K 


217.50 

0.230 

ro 


35.06 

0.156 

TABLE II: 

Ground state meson spectra, 

unit in MeV. 

States 

e 2 

A E 2 

PDG 

States 

e 2 

A E 2 

PDG 

7T 

142 

26 

139 

Vc 

2912 

5 

2980 

K 

492 

20 

496 

J/® 

3102 

4 

3097 

P 

826 

4 

775 

B° 

5259 

5 

5280 

CO 

780 

4 

783 

B* 

5301 

4 

5325 

K* 

974 

4 

892 

B° 

5377 

5 

5366 

<t> 

1112 

4 

1020 

B* s 

5430 

4 

5416 

D ± 

1867 

8 

1880 

B c 

6261 

7 

6277 

D* 

2002 

4 

2007 

Bt 

6357 

4 


Df 

1972 

9 

1968 

Vb 

9441 

8 

9391 

D* 

^ s 

2140 

4 

2112 

T(1 S) 

9546 

5 

9460 


The other details of the construction of the wave function 
can be found in our previous work E3- Subsequently, the 
converged numerical results can be obtained by solving 
the four-body Schrodinger equation 

= 0- ( 22 ) 

with the Rayleigh-Ritz variational principle. 

The energies E 4 ± A E\ of the charged states [cit][cd] 
with n 2S+1 Lj and J p under the assumptions of S = 
0,..., 2 and L = 0,..., 3 are systematically calculated and 
presented in Table III. The mass error of the states A E 4 
can be calculated just as A E 2 , they are around several 
MeV except for that of the state 1 1 Sq . The spin-orbit 
interactions are extremely weak, less than 2 MeV, there¬ 
fore the energies of excited states with the same L and 
S but different J are almost degenerate, see the ener¬ 
gies of the excited states with 1 5 I? 0 ) 1 5 -Di 1 5 -D 2 , 1 5 D 3 
and 1 5 H 4 in Table III, which is consistent with the con¬ 
clusion of the work [24]. Other spin-related interactions 
are stronger and bring about a larger energy difference 
than spin-orbital interactions, especially for the ground 
states with l 1 S'o, l 3 <Si and l 5 ^. The energy difference 
among excited states mainly comes from the kinetic en¬ 
ergy and confinement potential, which are proportional 
to the relative orbital excitation L. However, the rela¬ 
tive kinetic energy between two clusters [cu] and [cd] is 
inversely proportional to (X 2 ) while confinement poten¬ 
tial is proportional to (X 2 ) so that they compete each 
other to reach an optimum balance. 

The rms (r 2 ) 5 , (R- 2 ) 5 and (X 2 )a stand for the size 
of the diquark [cu ], the antidiquark [cd] and the distance 
between the two clusters, respectively, which are also cal¬ 
culated and listed in Table III. One can find that the 
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TABLE III: The energy E 4 + AE 4 and rms (r 2 ) 2 , (R 2 ) 2 and (X 2 ) 2 of charged tetraquark states [cm] [erf] with J p and n 2 S+ 1 Lj, 
unit of energy: MeV and unit of rms: fm. 


J p 

0+ 

0” 

0+ 

1+ 

1+ 

1" 

1“ 

1" 

1+ 

1+ 

n 2S+1 Lj 

VSo 

1 3 P 0 

1 5 A> 

l 3 Si 

2 3 S'i 

I’P, 

l 3 Pi 

l 5 Pi 

1 3 Di 

l 5 P>i 

E 4 i A_£/4 

3782 ± 12 

4097 ± 8 

4274 ± 7 

3858 ± 10 

3950 ± 10 

4075 ± 8 

4097 ± 8 

4153 ± 7 

4235 ± 7 

4273 ± 7 

(r 2 )i 

0.85 

0.96 

1.01 

0.90 

0.92 

0.94 

0.96 

1.00 

0.98 

1.01 

<r 2 >. 

0.85 

0.96 

1.01 

0.90 

0.92 

0.94 

0.96 

1.00 

0.98 

1.01 

<X 2 }^ 

0.42 

0.85 

1.12 

0.48 

0.66 

0.85 

0.85 

0.92 

1.10 

1.12 

J p 

1” 

2+ 

2" 

2“ 

2+ 

2+ 

2+ 

2" 

2" 

3" 

n 2S+1 Lj 

I s F x 

1 5 S 2 

1 3 P 2 

I s P 2 

l l D 2 

1 3 P> 2 

1 5 P> 2 

1 3 P 2 

1 5 F 2 

l 5 Ps 

_£/4 i A_£/4 

4387 ± 7 

4001 ± 7 

4096 ± 8 

4152 ± 7 

4212 ±8 

4235 ± 7 

4273 ± 7 

4354 ± 7 

4387 ± 7 

4150 ± 7 

to 

toll- 1 

1.02 

1.03 

0.96 

1.00 

0.95 

0.98 

1.01 

0.99 

1.02 

1.00 

(R 2 }2 

1.02 

1.03 

0.96 

1.00 

0.95 

0.98 

1.01 

0.99 

1.02 

1.00 

<X 2 }5 

1.30 

0.57 

0.85 

0.92 

1.09 

1.10 

1.12 

1.30 

1.30 

0.92 

j p 

3+ 

3+ 

3" 

3” 

3" 

4+ 

4~ 

4~ 

5" 


n 2S+1 Lj 

1 3 D 3 

1 5 T>3 

l'F 3 

1 3 F 3 

1 5 P 3 

1 5 T> 4 

1 3 P 4 

1 5 P 4 

1 5 P 5 


E 4 dz Ai ?4 

4234 ± 7 

4272 ± 7 

4332 ± 8 

4353 ± 7 

4386 ± 7 

4274 ± 7 

4353 ± 7 

4387 ± 7 

4387 ± 7 


(r 2 )i 

0.98 

1.01 

0.96 

0.99 

1.02 

1.01 

0.99 

1.02 

1.02 


<R 2 }5 

0.98 

1.01 

0.96 

0.99 

1.02 

1.01 

0.99 

1.02 

1.02 


<x 2 d 

1.10 

1.12 

1.27 

1.30 

1.30 

1.12 

1.30 

1.30 

1.30 



diquark [cm] and antidiquark [cd] share the same size in 
every Z+ state. The sizes of the diquark [cm] and an¬ 
tidiquark [cd] ((r 2 )2 and (R 2 ) 2 ) are mainly determined 
by the total spin S, the relative orbital excitation L of 
the states has a minor effect on them. However, the sizes 
do not vary largely with the total spin S, especially for 
higher orbital excited states. So the diquark [cm] and an- 
tidiquark [cd] are rather rigid against the rotation. For 
examples, the sizes of the two groups 1 1 S'o-1 3 5i- 1 5 52 
and changes gradually with the total spin 

S, 0.85-0.90-1.03 fm and 0.96-0.99-1.02 fm, respectively. 
And the sizes of the two groups l 1 £o-l 1 -Pi-l 1 -D 2 -l 1 -F 3 
and 1 3 5 i-1 3 Pi-1 3 Di- 1 3 F 2 vary slightly with relative or¬ 
bital excitation L, 0.85-0.94-0.95-0.96 fm and 0.90-0.96- 
0.98-0.99 fm, respectively. On the contrary, the distance 
between the diquark [cm] and antidiquark [cd] ((X 2 )s) 
changes remarkably with the relative orbital excitation 
L between the two clusters and is irrelevant to the total 
spin of the system, see the sizes of 1 3 S'i-1 3 Pi-1 3 Di- 1 3 F 2 
and l 1 5o-l 3 5’i-l 5 S' 2 _in Table III. The sizes of the diquark 
[cm], antidiquark [cd] and the distance between the two 
clusters are helpful to understand the changing tendency 
of energies of charged states Z + with quantum numbers 
S and L. 

In order to make clear the spatial configuration of 
charged states [cm] [erf], the distances in four states be¬ 
tween any two particles are given in Table IV. The ground 
state (l 1 ^ and 1+) of charged tetraquark [cm] [erf] pos¬ 
sesses a three-dimensional spatial configuration due to 
the competition of the confinement and the kinetic energy 
of the systems pa , which is similar to a rugby ball. The 
diquark [cm] and antidiqurk [cd] in the ground state have 
a large overlap because of the small (X 2 ) 2 , so the pic- 


TABLE IV: The average distances (r 2 j) 2 of the states [cu][cd] 
with 1 1 S'o, l 1 Pi, FPb, and PF 3 , r ij = n — rj, unit in fm. 


■ 2 S+ 1 Lj 

<r ? 2 >2 

( r § 4 ) 5 

( r 2 4 ) 5 

(r? 3 >^ 

( r ? 4 ) 5 

( r 23>^ 

<X 2 >2 

VS’o 

0.85 

0.85 

1.11 

0.46 

0.85 

0.85 

0.42 

l 3 Pi 

0.94 

0.94 

1.41 

0.87 

1.17 

1.17 

0.85 

VPb 

0.95 

0.95 

1.59 

1.11 

1.37 

1.37 

1.09 

1 1 p 3 

0.96 

0.96 

1.72 

1.28 

1.52 

1.52 

1.27 


ture of the diqurk or antiquark is not extremely distinct. 
However, all distances except for the sizes of the diquark 
and antidiquark ((r 2 2 )2 and (r| 4 ) 2 ) evidently augment 
with the increasing of the orbital angular momentum L 
in the excited states, see Table IV, which means that the 
picture of the diquark or antidiquark is more and more 
clear with the raising of the orbital angular momentum 
L. The spatial configuration of the excited states is still 
similar to a rugby ball, the higher orbital angular mo¬ 
mentum L, the more prolate of the shape of the excited 
states. The multibody color flux-tube basing on lattice 
QCD picture, a collective degree of freedom, plays an im¬ 
portant role in the formation of these charged tetraquark 
states, it should therefore be the dynamical mechanism 
of the tetraquark systems. 

Next, let us turn to argue the properties of the charged 
states Zj~ observed in experiments and their possible can¬ 
didates in the color flux-tube model, which are presented 
in Table IV. It can be seen from the Table IV that the 
spin and parity of the Z+( 3900) are still unclear up to 
now. The i?+(3900) may correspond to the same state 
as the Z+(3885) with 1+ [sj. The charged state [cu][cd] 
with 1 + and l 3 Si has a mass of 3858 ± 10 MeV in the 
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TABLE V: Z p states observed in experiments and their pos¬ 
sible candidates in the color flux-tube model. 




Experiment 



Model 


State 


Mass, MeV 

J p 

Mass, MeV 

J p 

n 28+1 Lj 

Z p (4050) 

m 

405111^4? 

? ? 

4075 ± 8 

1“ 

1 x Pi 

Z% (4250) 

yy 

4248^29-35° 

? ? 

4273 ± 7 

1+ 

1 5 Di 

Z p (3900) 

[2] 

3899.0t 3 ;«lt; 9 9 

? ? 

3858 ± 10 

1+ 

1 3 5i 

Z p (3885) 

[3J 

3883.9ijjl42 

1+ 

3858 ± 10 

1+ 

l 3 ^! 

Z p (3930) 

m 

3929^^2 

1+ 

3950 ± 10 

1+ 

2 3 S'i 

Zt (4025) 

[£] 

4026.3^;®!^ 

? ? 

4001 ± 7 

2+ 

l 5 S 2 

Z p (4020) 

m 

4022.9/°jt^ 

? ? 

4001 ± 7 

2+ 

1 5 52 

Z+( 4200) 

m 

4196i 3 ®/ 1 ^ 

? ? 

4235 ± 7 

1+ 

1 3 Di 

Zt (4475) 

m 

44751“!?® 

1+ 




o 

CO 

-f 

+ o 

N 

[9] 

44331?! 3 

1+ 





color flux-tube model, which is very close to those of the 
two charged states Z+(3885) and Z+(3900). It can not 
be excluded that the main component of Z+(3885) and 
Z+(3900) is the state [cu][cd] with 1 + and l 3 Si, which is 
supported by many theoretical work jl3]. The radial ex¬ 
cited state 2 3 Si has a mass of 3950± 10 MeV, which is ex¬ 
tremely close to that of Z+(3930). It is possible to iden¬ 
tify Z+ (3930) as the tetraquark state [cit] [cd] with 1 + and 
2 3 S'i. In other words, the Z+(3930) is the first radial ex¬ 
cited state of the Z+ (3900) in the color flux-tube model. 
The pair Z+(4020) and Z+(4025) show up with a sim¬ 
ilar mass (slightly above D*D* threshold). They might 
therefore be the same resonance, their spin and parity 
are unclear. QCD sum rule identified the Z+(4020) and 
Z+(4025) as a tetraquark state [cu] [cd] with 1 + [25j . 
the same approach also favored a tetraquark state but 
with different quantum numbers 2 + and 5 5 , 2 [26[. In 
our calculations, the nearest tetraquark state [cu][cd] to 
the Z+(4020) or Z+(4025) occupies quantum numbers 
2 + and l 5 ^. The tetraquark states [cu][cd] with 1“ and 
l 1 Pi and 1 + and l 5 Z)i have the energies of 4075 ±8 MeV 
and 4273±7 MeV, respectively, which are consistent with 
those of Z p (4050) and Z% (4250). So the two states may 
be assigned as the tetraquark states [cu] [cd] with 1“ and 
l 1 Pi and 1 + andl 5 Di, respectively, in the color flux-tube 
model. The newly observed Z+(4200) prefers 1 + , which 
can be described as the tetraquark state [cu][cd] with 1 + 
and l 3 Di in the color flux-tube model. The study of the 
three-point function sum rules on this state gives a sup¬ 
port to the tetraquark interpretation [24j. Of cause, it 
seems difficult to rule out other two possibilities of 2 + and 
1 1 I?2 and 2 + and 1 3 Z?2 in the model. The Z+(4430) is 
the first charged state, the J p of the state is determined 
unambiguously to be 1 + , the Zf( 4475) favors the spin- 
parity 1 + over other hypotheses J8|. Due to the heavy of 
the diquark and antidiquark, the energy of radial exci¬ 
tation between the diquark and antidiquark is too small 
to make the tetraquark state [cu] [cd] reach the energy 
above 4400 MeV. The internal excited states of the di¬ 
quark and/or antidiquark are needed to account for the 


heavy charged states, which is left for future. The al¬ 
ternative configuration for the two states may be meson- 
meson molecular states, which is suggested by several 
theoretical methods (l2j . 

From the above analysis and Table III, we can see 
that the most of the low energy theoretical states can 
be matched with the experimental ones. One of the ex¬ 
ception is the state with 0 + and RSo, which has a mass 
of 3780 ± 10 MeV. The experimental search of the r ) c -like 
charged state will give a crucial test of the present ap¬ 
proach. Our calculation also suggests that there are two 
negative parity states around 4100 MeV. More informa¬ 
tion on the the states around this energy is expected. 

The model assignments of the Z+ states are completed 
just in term of the proximity to the experimental masses, 
the more stringent check of the assignment is to study 
the decay properties of the states. These charged states 
should eventually decay into several color singlet mesons 
due to their high energy. In the course of the decay, the 
color flux-tube structure should break down first which 
leads to the collapses of the three-dimension spatial con¬ 
figuration, and then through the recombination of the 
color flux tubes the particles of decay products formed. 
The decay widths of the charged states [cu] [cd] are deter¬ 
mined by the transition probability of the breakdown and 
recombination of color flux tubes. The calculations are in 
progress. This decay mechanism is similar to compound 
nucleus decay and therefore should induce a resonance, 
which we called it as “color confined, multiquark reso¬ 
nance” state before [28[. 


IV. SUMMARY 

The charged tetraquark states [cw][cd] are system¬ 
atically investigated from the perspective of the color 
flux-tube model with a four-body confinement poten¬ 
tial. The investigation demonstrates that the charged 
charmoniumlike states Z+(3900) or Z+(3885), Z+(3930), 
Z+ (4020) or Z+(4025), Z+(4050), Z+(4250), and 
Z+ (4200) can be uniformly identified as tetraquark states 
[cu][cd] with the quantum numbers l 3 Si and 1 + , 2 3 Si 
and 1 + , 1 5 S 2 and 2 + , l 3 Pi and 1 _ , l 5 Z?i and 1 + , and 
l 3 Di and 1 + , respectively, in the color flux-tube model. 
The predicted lowest charged tetraquark state [cu] [cd] 
with 0+ and l 1 J So has a mass of 3780±10 MeV in the color 
flux-tube model. The model predictions would shed light 
on looking for possible charmoniumlike charged states in 
the future at the BESIII, LHCb and Belle-II. They fa¬ 
vor three-dimensional spatial structures which is simi¬ 
lar to a rugby ball,the higher orbital angular momentum 
L, the more prolate of the shape of the states. Those 
charged charmoniumlike states may be so-called “color 
confined, multiquark resonance”. However, the two heav¬ 
ier charged states Z+ (4430) and Z+ (4475) can not be de¬ 
scribed as tetraquark states [cu] [cd] in the color flux-tube 
model. 

The multibody color flux-tube is a collective degree of 
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freedom, which acts as a dynamical mechanism and plays 
an important role in the formation and decay of those 
compact states. Just as colorful organic world because of 
chemical bonds, multiquark hadron world should be var¬ 
ious due to the diversity of color flux-tube strucure. The 
well-defined the charged state Z+ (3900) and dibaryon d* 
resonance have been opening the gate of abundant mul¬ 
tiquark hadronic world. 
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